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SUMMARY 


A use is made of conformal mapping to determine equations that 
express the additional apparent mass of scalloped and/or clustered 
cylinder configurations in terms of the cross-sectional area of the 
configuration and the coefficients of a Fourier expansion of the loga- 
rithm of the radius vector of the configuration contour curve. Appli- 
cations are made to several configurations and comparisons with experi- 
mentally determined results are made for three of the configurations 
considered. The calculated and experimentally determined results are 
found to agree satisfactorily. 


INTRODUCTION 


A means for calculating the stability derivatives for a given body 
submerged in an incompressible fluid is through a use of apparent addi- 
tional masses or inertia coefficients that may be associated with the 
body as it undergoes various accelerated motions; see, for example, 
discussions on equations of motion and stability of bodies as presented 
in references 1 and 2. If the body under consideration is slender, 
certain of Its stability derivatives can be closely estimated by the 
well-known crossflow concept of Munk (ref. 3) and Jones (ref. 4) in 
which case the inertia coefficients of interest are those associated 
with a two-dimensional potential flow about the local cross sections 
of the body. (See refs. 5 to 7- ) Of current interest in connection 
with slender bodies are the inertia coefficients for sections of con- 
figurations resembling a cluster of circular cylinders, for example, the 
Saturn boosters, or perhaps cylindrical tanks with lengthwise corruga- 
tions that form a scalloped cross-sectional shape. These conf igurations 
will hereinafter be referred to as clustered and scalloped cylinders. 
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The inertia coefficients for a given configuration submerged in a 
perfect fluid medium can be directly related to either the kinetic 
energy imparted to the fluid or to the total forces that act on the 
configuration as it is accelerated in the fluid. This kinetic energy 
and/or these forces can be determined from a knowledge of the velocity 
potentials or, in the case of two-dimensional considerations, the com- 
plex potentials associated with the motions that the configuration 
undergoes. Moreover, in two-dimensional considerations the complex 
potentials, and hence the inertia coefficients, can be determined by 
means of techniques of conformal mapping which have been highly devel- 
oped in connection with airfoil theory; see, for example, references 8 
and 9* 

Expressions of the inertia coefficients for a variety of configu- 
rations, not including clustered or scalloped cylinders, however, can 
be found in references 1 to 7 . The purpose of this report is to make 
use of conformal mapping techniques to develop rather general expres- 
sions that will yield the inertia coefficients associated with trans- 
lations of either clustered or scalloped cylinders. These general 
expressions are then applied to an elliptic cylinder (scalloped cylinder 
with only two lobes), for a comparison with known results, and to two 
special configurations for comparisons with experimentally determined 
results. 


SYMBOLS 


A cross-sectional area 

A n , • * • constants associated with transformation of velocity 

potential in t^-plane (see, for example, eq. ( 13 )) 

t> semiaxes of an ellipse in £-plane 

^n constants associated with velocity potential for an arbi- 

trary cylinder (see eq. (l)) 

F a force per unit length acting at angle a 

F^ , F^ components of force per unit length parallel to 5- and 

rj-axes, respectively 

g acceleration due to gravity 

Z , L lengths 

Mq 


reference mass constant 
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M a apparent additional mass per unit length associated with 

accelerations in direction a 


in mass 

N number of circular cylinders involved in clustered 

conf igurat ion 

n integer 

R m radii 


R o> 

c 



constants 



, e 


polar coordinates in ^-plane 


^°, cp 

r 

t 

U, u, v 

V 

w 

x , y> S, *n 

z = x + iyl 

i = i + i 1 ]/ 

a 

6 = 0 - Cp 
7 \, U, 5 
A n C n , A n S n 
P 


polar coordinates in z -plane 

radius vector 

time 

velocities 

volume 

complex potential 
Cartesian coordinates 

complex variables 

angle measured counterclockwise from £- and x-axis 
parameters 

Fourier coefficients 

fluid density 

functions of 7\, N, Cp, 
and {6b) ) 


* 1 , x 2 


C 2 , and so forth (see eqs. ( 56 ) 


ANALYSIS 


A complex function of appropriate form is chosen to represent the 
complex potential associated with a cylinder of any cross-sectional 
shape (without circulation) and from a consideration of forces, and 
reactions thereto, acting on a section of this cylinder an expression 
for the apparent additional mass in terms of the cross-sectional area of 
the cylinder and certain constants appearing in the potential function 
is deduced. 


The Complex Potential Function and 
Apparent Additional Mass 

For convenience, the origin of coordinates Is chosen to he fixed at 
the centroid of a cross section of the cylinder and a fluid of constant 
density p is assumed to stream past the cylinder, normal to its axes, 
at some velocity U(t) from a direction that makes an angle a with 
the axis of reals (£-axis) as indicated in figure 1. It is further 
assumed that the cylinder is without circulation or rotation. Under the 
aforementioned assumptions the complex potential, w(^) will he of the 
form 


w(0 


-u(t) 



an + ih-j 

+ 


t 1^2 



+ 


^ + il3 n + 


( 1 ) 


where £ = | + irj is a complex variable, U(t) is the magnitude of the 
velocity at infinity and aj_,bi are constants that depend on the cross- 
sectional shape of the configuration and the direction of velocity a. 
The velocity components u and v parallel to the | and r\ axes, 
respectively, are given hy the relation 


dw „ 


u 


iv 


( 2 ) 


By a use of an extended theorem of Blasius (see, for example, sec- 
tions 9*52 and 9-53 of ref. 2) the components of force per unit length 
F| and F^ parallel to the |- and rj-axes, respectively, that act on a 

section of the cylinder may he expressed as 
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F| = -p^{2rtai - A cos a) (3) 

and 

F^ = -p^(2jtbi - A sin a) (4) 


where A is the cross-sectional area of the cylinder and a^_ and b-j_ 

are proportional to the real and imaginary parts, respectively, of the 
residue to the potential function. (See eq. (l) •) The force acting on 
the section of the cylinder in any direction in the plane of the section 
can, of course, he found hy appropriate projections of the components of 
force F| and F^. For example, the force F a acting through the 

centroid of the cylinder section and in the direction a of the stream 
is 


F 


a 


F| cos a + F Tj sin a =* 




cos a + bj_ 


sin a 



(5) 


It is of interest to note that a projection of forces normal to the 
direction a gives 


F/ Tf\ = F n cos a - Fj sin a = -2pjt^(hp cos a - ap sin a) (6) 

\ 2/ 1 ot x 


As will he seen at a later stage of this analysis, this normal component 
of force is generally zero for clustered and/or scalloped cylinders. For 
cases where it is not zero, it would cause the cylinder section to rotate 
unless a restraining force of equal magnitude were employed. 

Associated with the component of force acting in the direction a 
on the cylinder section is a reaction in the fluid, the net total of 
which corresponds to an opposing force. This opposing force may he 
expressed as a product of time rate of change in fluid velocity and an 
appropriate coefficient having the dimension of a mass per unit 

length 



-Fa 


( 7 ) 
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This mass term Mq, represents the additional apparent mass or inertia 
coefficient of interest in the present analysis. It may be expressed in 
terms of p, A, a x , U, and a by eliminating F a between equa- 

tions (5) and (7) • There is obtained 


Ma = 2jtp 



cos a + b-^ 


sin a 


) - 


A_ 

2jc 


(8) 


Equation ( 8 ) represents a key result in the present analysis. In this 
equation the cross-sectional area A for any given configuration can he 
considered as known. Thus the problem of determining becomes one 

of determining the residue ^ap + ibp) of the complex potential for the 

configuration under consideration. To this end the technique of con- 
formal mapping developed in connection with potential theory of arbitrary 
wing sections (refs. 8 and 9) are readily adaptable. 


Conformal Mapping as it Pertains to the Present Analysis 

As in potential theory for two-dimensional airfoils of arbitrary 
shape, the basic problem in the present analysis is to determine a 
conformal mapping function subject to restrictions pertaining to points 
at infinity that will transform a cylinder of arbitrary cross-sectional 
shape into a circular cylinder. The inverse of this mapping function 
will, of course, map the region outside the determined circle into the 
region outside the given closed area. Furthermore, this inverse mapping 
function will map a potential function for the circle into a corresponding 
potential function for the closed area. Hence, if one determines the 
mapping function that will transform a cylinder of given cross-sectional 
shape into a circular cylinder, the coefficients a n + ib n (eq. (l)) can 
be readily obtained from the inverse transformation and the noncirculatory 
potential function for a circular cylinder. Details of the method of 
determining the appropriate mapping functions are set forth in the fol- 
lowing section. 


The Mapping Function and Residue of the Complex Potential for a 
Clustered or Scalloped Cylinder Configuration 


Let a cross section of the given cylindrical configuration be 
assumed to be in the £-plane with its boundary defined by 
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5 = ( 9 ) 

where 9 = tan~l ^ is the polar angle, iy^( 0 ) = ^|2 + ^2 p s -^he 

radius vector to any poirrt on "the boundary and \|7 is a known function 
of 0. (See fig* 1.) The constant; R m is introduced simply to preserve 
dimensions and for definiteness can "be taken as the maximum radius of the 
configuration. Now let it be required to transform this configuration 
into a circle in the z-plane defined by 


z = R m e*0 +i * 


( 10 ) 


where cp is the polar angle and R^e^ (a constant) is the radius of 

the circle. The general transformation for this purpose is found in the 
developments of references 8 and 9 to be 


t ~ exp f(z) 


( 11 ) 


where f(z) is an infinite series of inverse powers of z. For clustered 
and/or scalloped cylinder configurations it is found convenient because of 
symmetry, to express this series as 


f(z) 


+ iB„ 


= “ n ■ 


n=l 


,nN 


( 12 ) 


where N is the number of cylinders in the cluster (or corrugations on a 
scalloped configuration) and A n and B n are constants, the determina- 
tion of which will be discussed along with the constant \|r 0 in the fol- 
lowing section. Note that at infinity the transformation defined by 
equations (ll) and (12) preserves the 1:1 correspondence of f and z 
'^1 


and that 


dz 


= 1 . 


By making use of equation (12) to expand the exponential term in 
equation (ll), the first few terms of the results may be written as 
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■where 


► n ^1 + ( A 2 + 5^2) + 1 ( B 2 + ^2) 

5 = Z 1 + zN + z2N 

( A 3 + c 3 + ^3) + ^ + D3 + ^3) 

+ 7 3N + * 


n 


J 


c 2 = A x 2 - B x 2 


D 2 = AiB x 


C J = a 1 a 2 " b 1 b 2 


Ex = A]_5 - 3A3B2/ 




D 3 = A 1 B 2 + B 1 A 2 
f 5 = b x 5 _ 3A 1 2 b 1 


J 


(13) 


(14) 


Inversion of equation (13) gives for the first few terms of z in terms 
of £ 


f Al + iBi A 2 + 2 + i[B 2 + (2N - 1)D 2 ] 

Z = ^\ 1 ' — S N " 5® 

A 5 + (3N - 1)C 5 + i(3H - 1) 2 Ex + i|B 5 + (3N - 1)D3 + UjN - 1) 2 P 3 

- 2 2 :=L + . . . > (15) 


,3N 


Nov the complex potential for a circular cylinder of radius 
out circulation is 



■with- 


v(z) 




(16) 


Substituting the expression for z given in equation (15) into equa- 
tion (l6) and expanding gives the potential for the cylinder configura- 
tion 
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These few terms of the complex potential are enough to determine the 
residue ap + ibp. Note the similarity of this potential function with 

that given in equation (l) and also note that the residue depends on the 
value of N; that is, 

for N = 1, 


ap + ibp = Rm e 


= R m 2 e 2ltr °e ia - 




for N = 


a l + ib l = R m 2e2 ^° eia " ( A 1 + iB l) e_ia 

= ^Rm 2e ^° " A ]) cos a ~ B 1 sin 

+ i ^R^e 2 ^ 0 + sin a - Bp cos a 1 


for N > 2, 


ap + ihp = R m 2 e 2 ^ < - ) e ia = R m 2 e 2 ^°( CO s a+isin a) 


( 20 ) 
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A clustered or scalloped cylinder configuration generally implies a value 
of N greater than unity. Thus, substituting values of ap and bp 
obtained from equations (19) and (20) into equation (8) gives: 


for N = 2, 


Mo, = 2rtp 


k 2 e% - A 


- Ap cos 2a - bp sin 2 a 


and for N > 2 , 


M a = 2 *P ba 2 ^ 0 ~ 


( 21 ) 


( 22 ) 


Observe that the expression for Mq, depends on the direction of motion a 

when N = 2 , but is independent of a when N > 2 . Also if the values 
of ap and bp, obtained from equation (20), are substituted into the 
expression for the component of force F rt normal to the direction of 

art— 

2 

motion (eq. (6)) it is found that this component of force vanishes when 
N > 2 . 


Determination of i|r 0 , Ap, and Bp 

By substituting the expression for t given in equation (9) and 
the expression for z given in equation (10) into the transformation 
given in equations (ll) and (12), there is obtained 


exp [t(e) - t 0 + i(9 - cp)J = exp 


00 

r A n 4* iBr^ / . 

— T^ cos “ 1 sin nJftp) 


1 M 


(23) 


By equating exponents on the two sides of this equation and then sepa- 
rating real and imaginary parts, there are obtained conjugate Fourier 
expressions for "the differences \|f - and 0 - cp, namely. 


t(e) - t 0 



1 



cos nNcp + B n 


sin nNcp) 


( 24 ) 
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0 - 9 =YZ 

i 






0 


nW 


(B n cos nNcp - A n sin riNcp) = e(cp) 


(25) 


It should be noted that through equation (25) 9 may be considered as a 

function of cp or, vice versa, cp may be considered as a function of 
0. It can now be recognized that and the coefficients A n and B n 

must have the following values (see, for example, ref. 8).: 



( 26 ) 


(27) 


( 28 ) 


It is known that two conjugate Fourier series, such as those in 
equations (24) and ( 25 ) , are related to one another through various 
forms of Poisson T s integral. A convenient form of this relation for the 
present purpose (see ref. 10) is 


€M ■ k 


■ 2rt 


i|/(9)cot 


9 ' - 9 


dtp' 


(29) 


This important relation, which may be looked upon as a functional 
equation relating >lrj^0(cp)j and c ( qp) is conducive to a rapidly con- 
verging iterative procedure for determining 1f{^(cp)] > and hence and 

the coefficients A n and B n , from equations (26), (27)* and (28). 

The scheme of iteration for ^^9(cp)J is simple in principle but 
involves some tedious expansions that are avoided in the examples in 
references 8 and 9 By recourse to graphical methods. The scheme is 
briefly outlined in the following paragraph and then applied to a general 
expression appropriate to clustered and scalloped cylinder configurations , 
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The Scheme of Iteration and Results for a Cylindrical 
Configuration of Arbitrary Shape 

It should be recalled that is assumed to be initially known as 

a function of 9 . As a first step in the iterative process then, the 
conjugate to t|t( 9) is determined by a use of equation (29) with 9 
replaced by 9. When the conjugate to t(9) is denoted by e-^9), the 

next step is to form the function t(9 + e-jJ and determine the conju- 
gate e 2 to this function. Then t(9 + e 2 ) is for m ed and its conju- 
gate is determined and so on. This process of iteration is logi- 

cally termed "successive conjugates" in references 8 and 9. It is 
demonstrated and emphasized in these references that only a few itera- 
tions or successive conjugates are usually necessary before ty(9 + e n J 

converges satisfactorily to i|fj^9(cp)J. This is borne out in the following 
applications . 

The logarithm of the radius vector R^^ ^ of cross-sectional 

shape of any clustered or scalloped cylinder configuration can be 
expressed as a trigonometric series of the form 



where Rq is a constant of linear dimensions and A, C n , and S n are 

nondimensional constants. The constant A, which may be a parameter 
peculiar to the configurations at hand or simply a chosen factor for the 
original coefficients of the trigonometric terms in the expansion of t, 
is factored out for convenience in the following developments. 

In view of the following relations (see, for example, ref. ll) 



cot a ~ — da = 0 


(31a) 
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2n 



cos ma cot 


a - 9 


2 


da = -sin m9 


and 


1 _ 

2rt 



sin ma cot 


SL-Z— da = cos m9 
2 


the conjugate Cp to ^(9) 


(eq. ( 301 )) is 



= El>P(Sn cos riN9 - C n sin nN9 


(311) 


( 31 c) 


(32) 


Following the procedure discussed in the foregoing paragraphs 
step is to form the function ^ ^9 + £pj . That is. 


C n cos nN^9 + e-^ + S n sin nN ^9 


/■D \ . 00 

*(e + Gp) = i°g(^) + XI 

v ' V^m / n=l 


the next 


+ 


(33) 


In order to obtain the conjugate e 2 to this expression the right-hand 
side may le expanded into a Taylor's series about = 0 and the 

results rearranged in powers of A. Terms of the expansion up to and 

p 

including A are 

t(9 + £p) = log ^ + ^(Cp 2 + Sp 2 ) + A(Cp cos N9 + Sp sin N9) 


+ A c 


(c 2 - |Cl 2 + |Si 2 )gos 2N9 
^S 2 - NSpCp^ sin 2N9 


( 34 ) 


Comparison of equation (j4) with equation (24) indicates that to a first 
approximation Ajp and are as follows: 


(35a) 


= 


R 


2/ 


log + ^f Cl 2 + Sl 




2 \ 


A 1 = A V R ra e 


t 


0 


N 


(35b) 


*0 


®1 “ A \®m e J 


N 


(35c) 


The expanded expression for the conjugate e 2 is found from equa- 
tion ( 3k) to he 

e 2 = A \S± cos N0 - Cp sin 0^ + A 2 ^S 2 - NSjC-^cos 2N0 

" ( C 2 " 7fl 2 + jjs^sin 2N0 + . . . 


(36) 


The next step is to form the function \Jr (0 + e 2 ) and then expand about 
e 2 = 0 and determine €j, and so forth. 

By extending the expansions of f(&+ e-jJ, \|^0 + } ... to 

include higher powers of A and proceeding as indicated in the fore- 
going^ all terms involving A to powers up to and including the sixth 
are found to have converged after three iterations; that is. 


t[e(cp)J = +(e + e 5 ) + o(a 7 ) 


(37) 


This relation indicates that, if A is small, the process converges very 
rapidly. Furthermore, the converged terms of ty(e + will generally 

be found to be sufficient for treating clustered or scalloped cylinder 
fv!^ would be , used for missi le boosters. Expressions 




A-^ and 


B 


*1 correct to the sixth power of A are found to 


involve at most six of the infinity of coefficients indicated in equa- 
lon (33), namely, C x , C 2 , Cj, S 1 , S 2 , and S*. The expressions 

are as follows: 
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Equations ( 38 ), (39) > and (^0) taken in conjunction with equations (21) 
and (22) constitute the main results of this analysis. The results are 
next applied to some specific conf igurations . 


APPLICATIONS AND COMPARISONS WITH EXPERIMENTAL RESULTS 


In order to conf inn the correctness of the foregoing analysis, appli 
cation of the results is first made to an elliptic cylinder for which an 
exact analytical expression for the apparent additional mass is known. 
Applications are then made to some cylinder shapes that would be of inter 
est in connection with clustered and scalloped cylinder configurations. 
Some of the results of these applications are then compared with experi- 
mentally determined results. 


Application to an Elliptical Cylinder 

The cross-sectional shape of an elliptic cylinder may be expressed 
In Cartesian coordinates as 



(a > b) (lH) 


In polar coordinates I = r cos 9, t] = r sin 9, the radius vector of the 
ellipse Is conveniently expressed as 


X (42) 

1 + - 2|i cos 29 


r = I^e ' = (1 + n)b- 

V 


■where 


Rm - a = 


1 + H 

1 - 


-\ 

-b 


l-i = 


a + b 


< 1 


(43) 


Taking the logarithm of each side of equation (42) gives 
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t(e) = log C 1 + jil^ _ 1 iog(l + p 2 - 2p cos 20) 

Km 2 


HO 


log — + u cos 20 + — cos 40 + . 
Rm 2 


u n 

+ — cos 2n0 + * 
n 


m 


By comparison with equation (50b) the various parameters are found to 
have the following values: 


*2.(i 

Em 


H = 2 

A = [JL 


C — i 


s n = 0 


- n) 






J 


(45) 


Substituting these quantities into equations ( 58 ), (59); and (40) gives: 

(46) 


t 0 = log(l - p) + p 2 + — + — 


= m^ 0 ) 




B-l = 0 




The cross-sectional area A of the cylinder is 


(47) 


A = Jtab 


(48) 


Substitution of the expressions for \|Tq, Ap, and A into equation (21) 
yields 
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M<x - 2;rtp/b 2 (l + p) 2 (l - n cos 2a)exp 2l|i 2 + — + — j - §£l> ( 49 ) 


A . ,6 


2 3 


It may be noted that for |i = 0, which implies that b = a, this 
expression reduces to 


Ma = *pa 2 (50) 


corresponding to the results for a circular cylinder of radius a. The 
exact expression for the apparent additional mass of an elliptic cylinder 

is (see ref. 7; page 239 for the equation for kinetic energy imparted to 
the fluid) 


M a = :rtp(a 2 sin 2 a + b 2 cos 2 aj ( 51 ) 

This expression may be written in terms of p as 



it may be seen that the argument of the exponential function in equa- 
tion (^9) is the same as the first few terms of this expansion. In fact 
for values of p less than about 0.5 (corresponding to values of eccen- 
tricity as large as 0 . 95 ), equation (49) gives essentially the same 
results as equation (51)* In order for the procedure to apply well to 
cylinder shapes with higher eccentricity, or values of p near unity, 
one would need to either extend the iteration discussed in a foregoing 
paragraph to include terms involving v to higher powers than the sixth 
or to employ an intermediary transformation, as discussed for example in 
references 8 and 9 in connection with airfoils, that would first trans- 
form the cylinder shape into a near circle. Such an intermediary trans- 
formation is not needed with the type of cylinder configurations con- 
sidered in the following paragraphs. 
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Application to a Family of Scalloped Configurations 

The equation of the radius vector of the periphery of the scalloped 
configuration considered may "be expressed as 

r = = R(l + C cos N0) (5*0 


■where R is a dime nsional constant, C is a nondimensional constant, and 
N ^ 5. A sketch of the configuration for N = 8 is shown in figure 2. 

In the Humeri cal examples that follow, the constant C is chosen so as 
to make the ratio of maximum radius vector to minimum radius vector a 
definite function of K. That is, for a given value of N, C is deter- 
mined from the relation 


1 + C = 
1 - C 



(55) 


where 


tan Aq 


/5 


l 


(56) 


Taking logarithms of equation (5^) and expanding gives 


\|r(0) = log — — h log(l+ C cos N9) 
Rm 


■p °° 

= log ^ + YZ *” c n 

% 1 


cos nN9 


(57) 


where 

A = 1 - A - eg 
c 

- s(-D n+1 

n S 

and 

R o _ c 

R m 2A(1 + C) 


(58) 

(59) 


( 60 ) 
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The cross-sectional area of the configuration under consideration is 
found to he 



' 2 rt 


0 


r 2 ae = 2^(2 + c 2 ) 


(61) 


By making use of equations ( 59 ) and ( 60 ) and noting in equation ( 57 ) that 
S n = 0 , the expression for i|Tq (eq. ( 58 )) pertinent to the configuration 

can he written as 


— log 7 — — r- + X -1 

0 2A(l + C) 1 


(62) 


where 


Xn = W 


2 A 2 + A 4 (l - 2N - 2N 2 ) + ^-(2 - 8N - 5N 2 + l 4 w 5 + mA) 


(63) 


The expression for Mq, (eq. (22)) thus becomes 




fi) 2 . 2 *! - (2 ♦ 02 ) 


(6b) 


It is convenient to refer the apparent additional mass of a given con- 
figuration to the apparent additional mass Mo of a circular cylinder. 

For this reference choose the cylinder with radius equal to the maximum 
radius of the configuration. In the present case 


Rm = R(1 + C) 


Hence 


and 


Mq = rtpR 2 (l + C ) 2 


M. 


‘a _ 1 


MO 2(1 + C ) 2 


2 2X / x 

-1 e 1 - (2 + C 2 ) 


'J 


( 65 ) 


( 66 ) 


(67) 
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This equation has teen evaluated for N = 3> • • • 8 and the results 

are shown plotted as discrete points in figure 3* As may he noted in 
this figure, the ratio of masses range from about 0 . 8 l for N = 3 to 
about 0.93 for N = 8 . About the same results are obtained for this case 
if the terms involving in (eq. ( 63 )) are ignored. Hence these 

results indicate almost complete convergence of the iterative process of 
determining \|Tq for this case. 


Clustered Cylinder Configuration 


The expression for x|r( © ) pertaining to a cluster of N circular 
cylinders can be obtained from a consideration of the equation of the 
part of the periphery pertaining to one of the N- cylinders in the 
cluster. For this purpose let the polar axis of the polar coordinate 
system be so chosen that the angular coordinate of the cylinder centers 

are — , — , . . . — ~ — rt . (See fig. 4.) The equation of the part of 
N N N 

the periphery pertaining to the cylinder with center at 9 = r = ^3 
may then be written as 


r = Rme^ 0 ) = Rjfcos^e - g) 


Note that for 
equation ( 68 ) 
interval 0 = 
c onf i gur at i on 



sin — ) . Taking logarithms of 
N / 


this case = Rj^l + 
and expanding the results into a Fourier series in the 
9 ^ £=£. gives the expression for i|f(9) pertinent to the 
under consideration. There is obtained 


*(9) 


-2 + C r 


ios sr 


r C n cos nN0 
nS 


( 69 ) 


where A is considered to be unity and the coefficients Cq and C n 
are determined as follows: 
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jt 



(see ref. 12 ) and 


it 



•where and denote Legendre functions (see ref. 13 ) . 

2 " T 

The expression for to for the present case can he written as 


= 



+ X 2 


where 


X 2 = | < Cp 2 + i( 8 C 2 2 + J+NCp 2 ^ - NpCp 4 ) + i. 288Cj 2 + ^NCpCpCj 

^ >j 

+ ^(cp^Cj - 3C! 2 C 2 2 ) - 56n 5 C 1 U C 2 + UN^J 


(72) 


( 73 ) 


The cross-sectional area of* the configuration under consideration is 
found to be 
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This expression has been evaluated for N = 3* 4 • • *8 and the results 
are plotted as discrete points in figure 5* The ratio of apparent addi- 
tional masses for this case may be noted in figure 5 to be slightly 
higher (ranging from about 0,86 for N = 3 to about 0,96 for N = 8) than 
results for the scalloped conf iguration in figure 3* Dropping the terms 
involving A^ in X 2 (that is, the terms in square brackets in equa- 
tion (73)) gives slightly lower values of the ratio M^Mq, especially 
for lower values of N. For example, for N = 3> the value of the ratio 
obtained when the terms involving A^ are not included is 0.840 whereas 
the value is O .865 when these terms are included. For N = 8, however, 
the two values are 0.957 and 0.964; thus there is almost complete con- 
vergence even when the A^ terms are not included in the value of X 2 . 


EXPERIMENTAL RESULTS AND COMPARISONS WITH THEORY 


Experimental Results 

Experimental values of the apparent additional mass of a scalloped 
cylinder conf iguration, a clustered cylinder configuration, and a single 
circular cylinder were determined by using water as a medium. The outer 
boundary of a cross section of the scalloped cylinder conf iguration 
(see fig. 2) is described by equation (54) with N = 8. The outer 
boundary of the clustered cylinder conf iguration (fig. 4) is described 
by equation (68) with N = 8. 

All the models were about 15.4 inches long and both the scalloped 
cylinder and clustered cylinder conf igurations had maximum radii of 
about 2.7 inches. The radius of the circular cylinder was 3-0 inches. 
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Each model was equipped with end plates and was suspended by 8 small 
cables arranged as shown in figure 6. The intended purpose of the end 
plates was to increase the effective length of the models and the cir- 
cular cylinder model was tested to obtain an estimate of the adequacy 
or effectiveness of the end plates. 

A cylindrical configuration suspended as indicated in figure 6 will 
behave as a simple pendulum with arm length l when each of the 8 cables 
is initially displaced through the same angle (3q (see view (b) in 

fig. 6) and released simultaneously. Hence, if the model is suspended 
in a fluid of density p and if all forces that contribute to drag and 

damping are neglected, the equation of motion may be expressed in terms 
of the angle p as 


(m + M a )ip = - g v(p c - P )p (76) 

where M denotes the mass of the model (cylinder plus end plate), 

denotes the total mass of fluid accelerated by the cylinder or the 
apparent additional mass, V denotes the volume of the model, p c 
denotes the mean density of the model, and g is the acceleration due 
to gravity. The appropriate solution to equation (7 6 ) is p = Pq cos cot 
from which one can obtain the following expression for M a 


Ma * - M (77) 

Thus, with a knowledge of g, V, p , p, l, M, and the measured 
frequency co of the pendulum one can determine Ma- In comparing 

experimental_re suits with calculated results, it was found convenient 
to express as a fraction 5 of the maximum radius of the 

cylinder as 


M, 


a 



(78) 


where L is the length of the conf iguration. Substituting this expres- 
sion into equation (77) and solving for 5 gives 
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(79) 


Comparisons of Experimental Results With Analytical Results 

As may he determined from the preceding section or from figures 5 
nnrl 5, the calculated results correspond to a value of 5 of about O.96 
for the scalloped configuration and about 0 .98 for the clustered cylinder 
configuration. The calculated value for the circular cylinder configura- 
tion is, of course, unity. 

Values of 5 determined experimentally and analytically for the 
three configurations tested are given in columns 8 and 9, respectively, 
of table I. The difference between the analytical and experimental 
values of 5 for the circular cylinder (l.OO and O.975) were used as a 
measure of the error in the experiment which was due to such things as 
the neglect of drag forces on both the configurations and the suspension 
cables, effect of wave motion, and possibly other factors in the equation 
of motion. Hence, the experimentally determined values of 5 for the 
scalloped cylinder configuration (0.938) and the clustered cylinder con- 
figuration (0.94l) were multiplied by the ratio & a nal/ 5 exp of the 
circular cylinder to obtain the adjusted experimental values which are 
presented in column 10 of table I. Comparison of the analytical and 
experimental values of 6 shows that the values obtained by analysis 
are a bit high but still in substantial agreement with measured values. 


CONCLUDING REMARKS 


It is concluded from the foregoing analysis and experimental evi- 
dence that the apparent additional mass of scalloped and/or clustered 
cylinder configurations can be accurately estimated from a knowledge of 
the cross-sectional area and the first few coefficients of a Fourier 
expansion of the logarithm of the radius vector of the configuration. 

It is remarked that the analytical procedure employed is not restricted 
to configurations for which equations of the periphery are known 
analytically because, with the use of harmonic analyses, the procedure 
can be readily adapted to graphically given configurations . 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May l4, 1962. 
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(a) Side view of suspended 
cylinder . 


(b) End view of suspended 
cylinder . 



(c) Cutaway view of setup. 

. - Sketch of experimental setup showing the clustered cylinder 
configuration with N = 8 and the suspension system. 
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